SOLAR CYCLE INDUCED VARIATIONS IN GONG p-MODE 
FREQUENCIES AND SPLITTINGS 
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ABSTRACT 

We have analysed the recently available GONG p-mode frequencies and splitting 
coefficients for a period of three and half years, including the rapidly rising phase of 
solar cycle 23. The analysis of mean frequency shift with different activity indices shows 
that the shift is equally correlated with both magnetic and radiative indices. During 
the onset of the new cycle 23, we notice that the change in 64 splitting coefficient is 
more prominent than the change in 62- We have estimated the solar rotation rate with 
varying depth and latitude. In the equatorial region, the rotation first increases with 
depth and then decreases, while an opposite behaviour is seen in the polar region. We 
also find a small but significant temporal variation in the rotation rate at high latitudes. 
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INTRODUCTION 



It is now well established that the frequency and splitting coefficients of the solar p-mode 
oscillations vary during the sun's magnetic cycle. The initial results were obtained from the study 
of the low degree modes arising from irradiance measurements (Woodard and Noyes 1985) and 
observations from the Sun as a star (Elsworth et al. 1990). Libbrecht & Woodard (1990) extended 
this study for higher £ modes ( 5 < £ < 140) and showed that the frequency shift increases with 
£ values. On the basis of observations made during the rising phase of cycle 22, Woodard et al. 
(1991) found a strong correlation between the frequency shift and the total surface magnetic field. 
The study of solar cycle 22 was extended to the initial phase of cycle 23 by several groups using the 
Global Oscillation Network Group (GONG) data and MDI data from the Solar and Heliospheric 
Observatory (SOHO). Anderson, Howe, & Komm (1998) analysed two 3-month equivalent data 
sets from GONG and MDI during the early phase of the solar cycle 23 and found evidence of small 
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shifts in frequency and 04 splitting coefficients. Using the GONG data, Bhatnagar, Jain, h Tripathy 
(1999) studied the variation in frequencies for a period of two years starting from 1995 August. 
They also presented evidence that magnetic and radiative activity indices correlate differently with 
frequency shifts. A similar result was earlier obtained by Bachmann & Brown (1993) for the 
declining phase of cycle 21 and ascending phase of cycle 22. Howe, Komm, & Hill (1999) have also 
analysed GONG frequencies over modes near 3 mHz and confirm the solar cycle changes in mean 
frequencies. In addition to frequency canges, they analysed the splitting coefficients up to a§ and 
reported that the even a coefficients show a larger correlation with the corresponding even Legendre 
component of the magnetic activity than with the average magnetic flux. Centroid frequencies and 
multiplet frequency splittings obtained from the MDI have also been studied by Dziembowski et 
al. (1998, 1999) who primarily concentrated on the solar cycle changes of the spherical structure of 
the Sun. 

Recently, Basu & Antia (1999) and Howe, Komm, & Hill (2000) found a small temporal varia- 
tion in GONG odd-order splitting coefficients upto 15th order. Using different inversion techniques, 
they found a systematic zonal flow migrating towards lower latitudes during the rising phase of the 
cycle. This flow can be associated with torsional oscillations as seen in BBSO data from 1986-1990 
(Woodard and Libbrecht 1993a). These zonal variations of the Sun's differential rotation were first 
noticed by Kosovichev & Schou (1997) using MDI data. They inferred subsurface flows from the 
rotational splittings of /-mode frequencies in the range of £ from 120 to 250 using 144 days time 
series. Toomre et al. (2000) have confirmed these zonal flows using both p- and /-mode splittings 
derived from MDI data. 

The changes in frequencies and splittings are believed to be caused by a combination of vari- 
ations in temperature, magnetic and velocity fields. Goldreich et al. (1991) proposed that the 
variations in solar eigen-frequencies are related to the perturbations in the surface magnetic fields. 
However, Kuhn (1998) argues that the frequency changes can not be due to the variation in surface 
magnetic field alone, because a reasonable explanation of the solar cycle acoustic changes must 
account for the apparent photometric cycle. 

In this paper, we analyse the frequencies and splitting coefficients obtained from the GONG 
project to extend previous works to all possible frequency and I ranges. We also infer the solar 
rotation rate as a function of depth and latitude using the analytical approach of Morrow (1988) 
involving the odd order splitting coefficients. Our results reveal a significant temporal variation in 
the rotation rate at high latitudes over the period of study. 

2. FREQUENCY SHIFT AND ACTIVITY INDICES 

The analysis presented here uses thirty three, 108-days frequency data sets (maximum available 
at the time of writing) from the GONG project, covering a period of more than 3 years from 1995 
May 7 to 1998 October 17. Each data set consists of 3 GONG months (GMs), where each GM 
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extends over 36 days and contains m-averaged p-mode multiplets. The mode frequencies were 
estimated using the standard GONG analysis (Hill et al. 1996) and obtained by fitting orthogonal 
polynomials as defined in Ritzwoller & Lavely (1991): 



where j s j(m) are the orthogonal polynomials for given value of degree I and c SyrL) i are the split- 
ting functions. The remaining symbols in Equation (1) have their usual meanings. Thus, this 
study differs from the earlier analysis of GONG data by Bhatnagar, Jain, & Tripathy (1999) and 
Howe, Komm, & Hill (1999), wherein the frequencies used were derived from Legendre polynomial 
expansion. 

The mean shift bv for a given £ and n is calculated from the relation 



where 8v n ^(t) is the change in the measured frequency for a given i and radial order n and a n ^ 
is the error in the observed frequency. The mean frequency shifts for GM 5-7 (1995 Sep. 23 - 
1996 Jan. 13) and GM 33-35 (1998 Jul. 2 - 1998 Oct. 17), which represents low and high activity 
periods respectively, are given in Table 1 where GM 1-3 is taken as the reference. It is noted that 
there are very few modes in the frequency range of 1500-1900 /iHz and 3900-4300 fiHz. It is further 
observed that the errors in 5v for these frequency ranges are larger as compared to other ranges. 

The variation of ^-averaged frequency shift for three different periods scaled with the inverse 
of the modal mass as a function of I is shown in Figure 1, where the frequency shift is calculated 
from the mean over all relevant data points. It is clear from the panels in Figure 1 that the shifts 
are constant for I between 40 to 80 in the frequency range of 1900-4300 /iHz whereas there is 
large scatter at both low and high degrees. The temporal evolution of the mean frequency shift for 
different I ranges scaled by the inverse of the modal mass is shown in Figure 2. Both these figures 
confirm that the frequency shift is a weak function of degree as was pointed out by Libbrecht 
& Woodard (1990). From GONG data, Howe, Komm, & Hill (1999) had also found no degree 
dependence for t approximately between 50 to 90 around 3mHz. 

In Figure 3, we show the variation of ^-averaged frequency shift as a function of frequency 
for three independent periods with reference to mean over all these periods. It is evident that 
the frequency shift has a strong dependence on the frequency which increases with increasing 
frequencies. This is also consistent with the result of Libbrecht h Woodard (1990) and Anguera 
Gubau et al. (1992). The temporal evolution of frequency shifts over the entire period of three and 
half years for different frequency ranges is shown in Figure 4. It is noted that the frequency shifts 
for all frequency ranges are small upto GM 24-26 and then increases rapidly. This behaviour of 
the shifts can be understood in terms of the rapid increase of solar activity represented in the same 
figure by means of the 10.7 cm radio flux. It is further noticed that the frequency shifts for all 
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frequency ranges seem to converge at GM 23-25. This particular aspect was investigated in detail 
using different combination of GMs as reference months. We conclude that this converging feature 
depends on reference frequencies and has no physical significance. 

The mean shift is correlated with different solar activity indices by computing a mean value for 
each index over the same 108 days interval over which the frequencies were determined. The activity 
indices considered are: Ri, the International sunspot number obtained from the Solar Geophysical 
Data (SGD); KPMI, Kitt Peak Magnetic Index from Kitt peak full disk magnetograms; SMMF, 
Stanford Mean Magnetic Field from SGD; MPSI, Magnetic Plage Strength Index from Mount 
Wilson magnetograms; FI, total flare index from SGD and T. Atag (1999, private communication); 
He I, equivalent width of He I 10830 A line, averaged over the whole disk from Kitt peak; Mg II, 
core-to-wing ratio of Mg II line at 2800 A from L. Floyd (1999, private communication); F±q, 
integrated radio flux at 10.7 cm from SGD. A linear least square fitting is carried out between 
mean activity indices and frequency shifts. The fitting parameters, \ 2 values and both Pearson's 
(r p ) and Spearman's (r s ) correlation coefficients are summarised in Table 2. It is observed that 
both the magnetic and radiative indices are equally correlated with the frequency shifts in contrast 
to the earlier results of Bachmann & Brown (1993) and Bhatnagar, Jain, & Tripathy (1999). This 
aspect needs to be investigated in detail with the availability of more high precision data sets in 
future. The correlation coefficients for different frequency ranges are given in Table 3. It is observed 
that the correlation has an increasing trend with frequency ranges and are marginally different at 
higher frequencies. 



3. SPLITTING COEFFICIENTS 

It is known that the solar differential rotation and other symmetry breaking factors like mag- 
netic field can lift the degeneracy of the solar acoustic modes and split the eigen frequencies. 
Individual mode splittings can be represented by Legendre polynomial expansion in (m/L) 

m 

+ L a s , n /Ps {m/ L) , (3) 

s=l 

where the expansion coefficients a SjTl ^ are known as the splitting coefficients and I? = £(£ + 1). 
Alternatively, the splitting data can be expressed in terms of the splitting functions, c SjTl ^, as given 
in Equation (1). These functions are generally converted to a new set of splitting coefficients 6's by 
using the relation (Ritzwoller & Lavely 1991) 

_ w v 

where 

*; = (-D<») (; { jj, (5) 
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and 



{2i-s)\ 



.1/2 



(2^ + s + l)! 



(6) 



The odd-order coefficients measure the solar rotation while the even-order coefficients probe the 
symmetry about the equator. The nonzero values of these even coefficients reflect the pole-equator 
asymmetries in the solar structure. It may be noted that we have retained terms up to s = 5 
because the coefficients and their errors are of similar magnitude for s > 5 in the data sets used in 
this study. 

The temporal variation of 62 and 64 coefficients is shown in Figure 5. We find that 62 has strong 
correlation with activity, whereas 64 is anticorrelated. The results of linear regression analysis 
between these coefficients and activity indices are given in Table 4. It may be noted that 62 has a 
better quantitative correlation than 64. We observe that during the onset of cycle 23, the change 
in 64 is more prominent than the change in 62. Howe, Komm, & Hill (1999) have also studied the 
temporal variation of the 62 and 64 coefficients at 3 111Hz and reported a phase shift of one year 
between them. The odd order coefficients, 61, 63 and 65, calculated using the expression (4) are 
used to obtain information on the variation of the Sun's rotation rate with depth and latitude. 



The solar rotation rate using the helioseismic data are generally obtained through two different 
methods: in forward approach the frequency splittings are computed for a chosen solar rotation 
model and then compared with the observed splittings. In the inverse method the measured fre- 
quency splittings are used directly to produce a single function for the angular velocity. In this 
study, we use the analytical method of Morrow (1988) where the appropriate combination of odd 
order splitting coefficients reflects the depth variation of angular velocity at chosen co-latitude (p, 

o 

(90 — 8, where 6 is the latitude). 



3.1. 



Solar Rotation Rate 



2 




(7) 



i=0 



where 



di = 1, 

c?3 = [1 — 5cos 2 (f)], 

d 5 = [1 - Ucos 2 (j) + 21cos 4 <p] 



(8) 
(9) 
(10) 



The corresponding rotation rate at equator is given by 



n nl (90° ) = af + af + af, 



(11) 



and at pole 



(0°) = af -4af + 8af. 



(12) 
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Wilson, Burtonclay & Li (1997) have recently modified expressions (7-12) by including higher 
order splitting coefficients up to 07. The addition of further terms to these expressions seems to be 
of limited utility due to the larger relative errors of higher order coefficients. We compute a's from 
c's using Equations (61-63) of Ritzwoller & Lavely (1991). In Figure 6, we show the variation of 
ai, (Z3 and CI5 with degree. The non-zero values of these coefficients indicate that the solar rotation 
rate changes with depth and latitude. 

In Figure 7, we show the sun's rotation rate at different latitudes as a function of depth 
approximated by v/L which corresponds to the radius of lower turning point r t as defined by the 
relation (Christensen-Dalsgaard & Berthomieu 1991) 



where c is the sound speed. A higher value oiv/L denotes a smaller value of r t and hence a greater 
depth. As a result the low degree modes are sensitive to the rotation from surface to the core, while 

o 

the high degree modes probe solar rotation near the surface. From Figure 7, we find that for 9 = 

o 

and 30 , the rotation rate first increases and then decreases as depth increases. On the other hand, 

o o 

the rotation rate for 9 = 45 and 60 remains constant below the base of the convection zone (v/L 

o 

Ri 70), beyond which it starts increasing. We further notice that for 9 = 60 , the increase in internal 

o 

rotation rate above v/L = 70, is steeper than at 9 = 45 . This suggests that, as the depth increases, 
the equatorial rotation rate decreases while the polar rotation rate increases. Wilson, Burtonclay 
& Li (1997) have analysed frequency splittings derived from first four months of GONG data and 
showed that the greatest changes in rotation occurs in the region below the base of the convection 
zone. They also found that the data does not support models exhibiting a discontinuous shear 
between the convection zone and a uniformly rotating radiative region. Thompson et al. (1996) 
have also estimated the solar rotation rate using the first few month's of GONG data. Our results 
derived from analytical expressions are in gross agreement with earlier studies. 

The observed rotation rate is conventionally expressed in terms of even powers of cos</>: 



where (j) as defined earlier is the solar co-latitude. The surface constants A, B, and C are related 
to splitting coefficients by the relations 



tt(r, <t>) = A + B cos 2 (f> + C cos 4 (f>, 



(14) 




(15) 



(16) 



(17) 



Various workers have calculated these surface coefficients using different data sets. Brown et 
al. (1989) found A = 462.8 nHz, B = - 56.7 nHz and C = - 75.9 nHz for r > 0.723 R Q using 
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Call K intensity data taken from South pole. Based on 100-day observations made at BBSO for 
£ between 10 to 60, Libbrecht (1989) found the best fit with A = 461 nHz, B = - 60.5 nHz and 
C = — 75.4 nHz. We have derived these constants from the GONG data for the 20 < I < 100 
and obtained the best fit with A= 459.54 nHz, B = - 61.90 nHz and C = - 70.20 nHz. Using 
the odd-order splitting coefficients of GONG data derived from Legendre polynomial expansion, 
we obtain, A = 459.97 nHz, B = — 61.33 nHz, and C = - 70.93 nHz. The values obtained for A, 
B and C using GONG data are in close agreement with the earlier values. 

In Figure 8, we plot the average surface rotation rate as a function of latitude using the derived 
coefficients. In the same plot, we have also shown the results from inversion techniques (Antia, 
private communication) and Doppler surface measurements (Snodgrass 1984). It is clear that the 
rotation rate changes significantly in mid latitude while the change in rotation rate near pole and 
equator is small. We find that the rotation rate derived from the GONG data agrees well with 
other results. However, the inverted rotation profile departs from other rotation rates beyond the 

o 

latitude of 70 , probably due to the resolution limitation in inversion techniques. 

Woodard and Libbrecht (1993b) studied the time variation of equatorial rotation rate for differ- 
ent £ ranges and suggested a small variation from year to year. To investigate temporal variations, 
Figure 9 shows the residual surface rate at four different latitudes, obtained after subtracting the 
average angular velocity shown in Figure 8 by long dash line. It is evident that at higher latitudes, 
the residual rotation rate, commonly known as zonal flows, is time dependent. This has a magni- 

o 

tude of approximately 2 nHz at a latitude of 60 . Recently Basu & Antia (1999), Howe, Komm, 
& Hill (2000) using GONG data, and Toomre et al. (2000) using MDI data, have also reported a 
small but significant time variation in the rotation rate. Our results derived from the analytical 
method are consistent with these inversion studies. 

4. CONCLUSION 

The mean frequency shift varies over the solar cycle and is correlated similarly with both the 
magnetic and radiative indices. 

The &2 splitting coefficient is linearly correlated with the activity indices, while 64 is anti- 
corr elated for all degree and frequency ranges. 

We detect a small but significant variation in the rotation rate derived from the linear combi- 
nation of odd order coefficients over a period of three and half years. 

We thank the anonymous referee for his critical remarks which impreved the manuscript. We 
also thank T. Atag, L. Floyd, and R.K. Ulrich for supplying us the Flare index, Mg II, and MPSI 
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obtained by the Global Oscillation Network Group project, managed by the National Solar Obser- 
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Fig. 1. — The binned frequency shifts for three different GONG months, scaled by the inverse of 
the modal mass, as a function of degree i. The shifts have been calculated with respect to the 
mean of all three periods. Error bars represent mean error in shift. 

Fig. 2. — The temporal evolution of the frequency shifts scaled by the inverse of the modal mass for 
three different t ranges: 20-40 (traingles), 40-80 (stars) and 80-100 (diamonds). The corresponding 
period is given in the upper x-axis. 

Fig. 3. — The binned frequency shifts for three independent periods: GONG months 8-10 (squares), 
23-25 (diamonds) and 33-35 (triangles) as a function of frequency. The shift has been calculated 
with reference to the mean of all three periods. Error bars represent mean error in shift. 

Fig. 4. — The temporal evolution of frequency shifts for different frequency ranges: 1900-2300 fiHz 
(plus), 2300-2700 fiRz (stars), 2700-3100 fiRz (squares), 3100-3500 ^Hz (dimonds), and 3500- 
3900 fiHz (triangles). The solid line without symbols represents the total mean shift in the range 
of 1900-3900 fiKz. The shifts have been calculated with respect to the mean of all 33 data sets. 
The solid line with crosses represents the scaled 10.7 cm radio flux. 

Fig. 5. — The temporal evolution of 62 (solid line) and 64 (dashed line) coefficients derived from 
the splitting functions c s . The solid line with diamonds represents the scaled 10.7 cm radio flux. 

Fig. 6. — The variation of splitting coefficients with degree I for GM 8-10. The solid line represents 
the average values of the coefficients for t = 20-100. The la errors are shown by dotted line. 

Fig. 7. — The solar rotation rate as a function of vjL, a proxy for the turning point radius r t 

OOO o 

(Equation 13), for five different latitudes: 9=0 (plus), #=30 (stars), 0=45 (diamonds), 9=60 

o 

(triangles), and 9=90 (squares) are shown. 

Fig. 8. — The variation of average solar rotation rate calculated from three term fits in even 
power of cos0 are shown as a function of latitude. The solid line represents the rotation rate from 
Snodgrass (1984), dotted and dash-dot-dot-dot lines are for BBSO data sets, as derived by Brown 
et al. (1989) and Libbrecht (1989) respectively. Other lines represent rotation rates calculated from 
GONG data; long dash line for analytical results using orthogonal polynomials, and short dash for 
inversion results (Antia, private communication). 

Fig. 9. — The temporal evolution of residual rotation rate near the surface at four different latitudes. 
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Table 1. NUMBER OF COMMON MODES AND FREQUENCY SHIFTS FOR I BETWEEN 

20 TO 100 AND DIFFERENT v RANGES 



Frequency range Number Frequency shifts 



(//Hz) of modes (nHz) 











GM 5-7 


GM 33-35 


1500 < 


V 


< 1900 


6 


11.4 ± 5.51 


32.2 ± 6.12 


1900 < 


V 


< 2300 


71 


-04.4 ± 1.71 


40.9 ± 1.86 


2300 < 


V 


< 2700 


92 


-05.8 ± 1.64 


96.5 ± 1.82 


2700 < 


V 


< 3100 


66 


-16.7 ± 1.71 


173.7 ± 2.06 


3100 < 


V 


< 3500 


84 


-13.4 ± 1.99 


270.4 ± 2.39 


3500 < 


V 


< 3900 


110 


-29.3 ± 3.63 


338.3 ± 4.05 


3900 < 


V 


< 4300 


31 


-34.5 ± 16.4 


377.5 ± 17.6 


1900 < 


V 


< 3900 


423 


-10.4 ± 0.08 


139.1 ± 0.10 
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Table 2. FITTING AND CORRELATION STATISTICS IN THE FREQUENCY RANGE OF 

1900-3900 fiHz AND t BETWEEN 20 TO 100 



Activity index 


slope 


intercept 


x 2 


r p 


r s 




(nHz per activity unit) 


(nHz) 








Ri 


2.154 ± 0.07 


-52.802 ± 0.24 


14.0 


0.98 


0.93 


KPMI 


25.04 ± 0.90 


-205.67 ± 0.74 


26.6 


0.96 


0.82 


Fio 


2.758 ± 0.09 


-230.84 ± 0.83 


10.3 


0.99 


0.95 


He I 


8.416 ± 0.30 


-391.20 ± 1.41 


31.5 


0.95 


0.92 


SMMF 


7.283 ± 0.36 


-85.356 ± 0.44 


397 


0.87 


0.69 


MPSI 


136.6 ± 4.87 


-48.772 ± 0.23 


9.20 


0.99 


0.92 


FI 


27.69 ± 1.58 


-39.642 ± 0.22 


50.4 


0.93 


0.88 


Mg II 


19580 ± 697 


-5001.6 ± 17.8 


7.30 


0.99 


0.93 
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Table 3. CORRELATION STATISTICS FOR DIFFERENT FREQUENCY RANGES* 



Activity index 1900-2300 2300-2700 2700-3100 3100-3500 3500-3900 

p <s <p 's >p >s >p >s <p 's 



Ri 0.92 0.86 0.97 0.88 0.98 0.92 0.99 0.95 0.98 0.93 

KPMI 0.89 0.76 0.94 0.74 0.96 0.81 0.98 0.89 0.98 0.88 

Fio 0.91 0.89 0.97 0.89 0.98 0.95 0.99 0.97 0.99 0.96 

He I 0.84 0.86 0.93 0.84 0.95 0.90 0.97 0.97 0.97 0.97 

SMMF 0.70 0.69 0.77 0.68 0.80 0.67 0.82 0.68 0.82 0.72 

MPSI 0.93 0.87 0.98 0.86 0.99 0.92 0.99 0.95 0.99 0.96 

FI 0.84 0.84 0.91 0.85 0.92 0.88 0.94 0.90 0.94 0.87 

Mgll 0.92 0.86 0.98 0.86 0.99 0.93 0.99 0.97 1.00 0.96 



* Frequencies are in [iKz 
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Table 4. FITTING AND CORRELATION STATISTICS FOR b 2 AND 6 4 COEFFICIENTS 



Activity 


slope 


intercept 


x 2 






index 


(nHz per activity unit) 


(nHz) 










Ri 


10.26 ± 3.09 


-221.48 ± 93.60 


1.0 


0.92 


0.80 


KPMI 


122.28 ± 36.19 


-973.70 ± 298.10 


0.6 


0.95 


0.82 


^10 


13.22 ± 3.94 


-1076.61 ± 330.48 


0.8 


0.93 


0.81 


He I 


40.96 ± 12.21 


-1873.15 ± 565.20 


0.8 


0.94 


0.88 


SMMF 


34.80 ± 14.60 


-378.15 ± 174.85 


6.3 


0.88 


0.75 


MPSI 


656.86 ± 195.16 
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